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SHORELINE AS A CONTROLLING FACTOR 


IN COMMERCIAL SHRIMP PRODUCTION 


By Kenneth H. Fallen 
National Space Technology Laboratories* 


SUMMARY 


The coastal zone of the United States is an area subject to tremendous pressures, 
as population centers expand and the impacts of residential, agricultural, industrial, 
and commercial factors propagate through coastal wetlands, bays, and estuaries. 
Federal and, in many instances, state and local legislation requires management and 
monitoring of coastal resources, including the most basic life forms and processes. 

In response to these requirements, efforts are being made by many to develop a de- 
tailed understanding of basic coastal processes and the influence of man's activities on 
these processes, and to develop techniques for monitoring them. The research effort 
described in this work was conducted to develop an understanding of the importance of 
a single process in the overall system and to demonstrate a technique by which a con- 
trolling factor in that process can be monitored synoptically using satellite data. 

Using the systems ecology approach, it is possible to trace the flow of energy 
from the marsh ecosystem to the estuary ecosystem, and to relate secondary produc- 
tion in the estuary to this influx of energy. The influence of shoreline as a limiting 
factor on the flow of energy-carrying nutrients from the marsh to the marine ecosys- 
tem, as reflected by the commercial harvest of shrimp in Louisiana bays and estu- 
aries, was analyzed. Data acquired by the Landsat multispectral scanner were com- 
puter-processed to develop statistics relating to detritus production in the marshlands 
and the length of the marsh-water interface. These statistics were found to have a 
highly significant correlation with the commercial shrimp harvest. The data were 
used to develop a mathematical model based on detritus production and export to pre- 
dict the long-term average commercial shrimp harvest for nine segments of the Lou- 
isiana coast. Detritus production was estimated to be proportional to the area having 
more than 5.0 kilometers of shoreline per square kilometer of area (5.0 km shore- 
line/km 2 ) and export to the marine ecosystem was modeled as the product of shoreline 
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length and mean tidal range. The result was an excellent agreement between reported 
and predicted long-term average harvest, with the root mean square deviation between 
the reported and predicted values being 4.36 kg/ha over a range of 0. 29 to 45. 16 kg/ha. 
The analysis thus indicates that the production of detritus on the marshlands and its ex- 
port, as regulated by the tidal flow across the shoreline, are conti’olling factors in the 
production of shrimp in the Louisiana bays and estuaries. 

With further research, it should be possible to extend the analysis to provide an 
important tool for coastal zone management. Remote sensing can be used to monitor 
marshlands and routinely assess biogeographic al factors. Trends of changes taking 
place in the marsh, whether natural or anthropogenic, and proposed modifications to 
the marsh could be analyzed with an ecosystem model similar to the one developed in 
the present work to forecast possible changes in long-term average shrimp production. 


INTRODUCTION 


The coastal zone of the United States is an area subject to tremendous pressures, 
as population centers expand and the impacts of residential, agricultural, industrial, 
and commercial factors propagate through coastal wetlands, bays, and estuaries. Fed- 
eral and, in many instances, state and local legislation requires management and mon- 
itoring of coastal resources, including the most basic life forms and pi'ocesses. In 
response to these requirements, efforts are being made by many to develop a detailed 
understanding of basic coastal processes and the influence of man's activities on these 
processes, and to develop techniques for monitoring them. The research effort de- 
scribed in this paper was conducted to develop an understanding of the importance of a 
single process in the overall system and to demonstrate a technique by which a con- 
trolling factor in that process can be monitored synoptically using satellite data. 

The advent of systems ecology has made possible the analysis of the various com- 
ponents of natural ecosystems. While adequate data are seldom available for complete 
mathematical treatment of a model ecosystem, it is still possible to assess the sig- 
nificance of individual elements. To model an ecosystem, one identifies the subsystems 
which can be separated as discrete entities and the processes and paths of energy flow 
relating them. As a first-level analysis, the marsh-estuary system may be divided 
into two major subsystems, terrestrial and aquatic, linked by the flow of organic and 
mineral nutrients carrying chemical potential energy. The principal transport mech- 
anism linking the terrestrial and aquatic systems is the flow of water across the shore- 
line under the influence of tidal fluctuation and rainfall. It is generally recognized that 
an important factor contributing to the tremendous productivity of salt marsh estuaries 
is the interaction between the marsh and the water (e.g. , Schelske and Odum, 1961; 
Teal, 1962; Day et al. , 1972). This paper presents an analysis of the influence of 
shoreline as a limiting factor on the flow of energy-carrying nutrients from the marsh 
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to the bays and estuaries as reflected by the commercial shrimp harvest in the Louisi- 
ana estuaries and bays. This analysis is based on data derived from published sta- 
tistics relating to the marsh-estuary biology and from computer analysis of satellite 
mappings of the Louisiana coast. 

Mapping of coastal wetlands is a very difficult problem. In addition to the tremen- 
dous difficulty of performing field surveys in the wetlands, these areas are subject to 
constant change. Maps prepared from data acquired during 1950-1960 show significant 
deviations from current aerial photography. The use of multispectral scanner (MSS) 
data from the Landsat satellite ameliorates the problem by providing the capability to 
routinely monitor the wetlands, to update existing maps, or to generate original maps 
based on identifiable control points located on existing maps. From computer analysis 
of the data acquired by the Landsat MSS, thematic maps showing land and water, vari- 
ous species of vegetation, and residential or industrial development can be produced. 
Computer processing of the Landsat MSS data, available directly in computer-compat- 
ible form, makes feasible the routine monitoring of extensive areas, such as the entire 
coast of a state. 

The research upon which this report is based required the mapping of nearly the 
entire coast of Louisiana, a task that would have been impossible by any conventional 
techniques within the constraints of reasonable funding. Analysis based on existing 
maps would have been questionable due to significant changes that have taken place in 
the coastal wetlands since the maps were produced and significant errors in the initial 
mapping. The Landsat MSS data provide the opportunity to develop geographic param- 
eters over very large areas, with good accuracy, at a reasonable cost. 


THEORY 


Ecosystem models for marsh-estuary environments have been developed by vari- 
ous researchers, including Carter et al. (1973), Teal (1962), and Day et al. (1973). 

Each of these models emphasizes the importance of the link between terrestrial and 
aquatic subsystems. Figure 1 is the simplified energy flow diagram of the marsh- 
estuary system in Barataria Bay published by Day et al. (1973). Tide, water level, and 
rainfall are important forcing functions which drive the flow of inorganic nutrients, salt, 
and detritus between the subsystems. Nutrients carried by river waters find their way 
into the estuary under the influence of tidal action and fertilize the marsh, whereas 
under the same influence, detrital materia 1 is washed from the marsh into the estuary 
and eventually the Gulf of Mexico. Human involvement occurs with the harvest of estu- 
arine fauna and discharge of waste material. This general form developed for Barataria 
Bay is applicable to the entire Louisiana coast, the study area for the subject analysis. 

Figure 2 is a detailed schematic of the portion of the ecosystem centered at the 
land/water interface, developed to show the energy flow leading to the only consumer 
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Figure 2. - Energy flow diagram leading to in-shore shri 








studied in tne subject analysis, penaeid shrimp. The marine subsystem is further sub- 
divided into the bay-estuary subsystem, consisting of semi-enclosed water bodies and 
interconnecting bayous and channels; the coastal subsystem, including the open waters 
along the coast, outside the bays and sounds; and the deep Gulf. These shrimp, to- 
gether with other species including amphipods, mysids, oslracods, planktonic cope- 
pods, crabs, filter-feeding bivalves, and a few species of fishes, arc detritus con- 
sumers; they derive a significant amount of their nourishment from the ingestion of 
vascular plant detritus together with small quantities of live algae (Odum et al. , 1972). 
Most of the detritus available to the shrimp finds its way into the bays and estuaries 
from the marsh, as tidal flow and rainfall wash dead plant material from the marsh 
subsystem into the bay-estuary subsystem. For one bay in Louisiana, Day ct al. (1973) 
estimated that as much as 70 percent of the total organic production available in the 
water was detritus from marsh grasses. The remainder is produced by plankton and 
benthic vegetation directly in the bay subsystem. Some detritus is riverborne, and is 
carried into the bays from the coastal waters by tidal action and currents. The river 
water is rich in inorganic nutrients, which are also carried into the bays by the tides. 
There is evidence that potassium, magnesium, and phosphate from these waters fer- 
tilize the marshland, whereas nitrates appear to be leached from the l;md by rain and 
tidal flow (Palmisano, 1970). River discharge also regulates the salinity of the coast- 
al and bay waters. As indicated by figure 2, the model is based on the hypothesis that 
the production of shrimp (a self-maintaining consumer) represents a direct flow of 
energy from detritus (an active energy storage factor) when salinity and temperature 
are in the proper range. If this hypothesis is correct, shrimp production should be 
related to the production and transport of detritus from the marsh into the bay. 

The transport of detritus into the bay from the marsh is controlled by two work 
gates, labeled 1 and 2 in figure 2, which operate under the influence of tidal action and 
the runoff of rainfall, respectively. The first gate is bidirectional, whereas the second 
permits rain-induced flow only from the marsh to the bay. Some detritus is carried 
into the coastal waters by the rivers, and consequently is transported into the bays by 
incoming tides. The transport of this terrigenous detritus is controlled by work gate (>. 
Similar work gates control the flow of inorganic nutrients and salt. 

The form of the mathematical function describing these work gates can be defined 
to a first approximation by synthesizing the significant factors influencing the trans- 
port mechanisms. The first factor to consider is the "conductivity" of the interface 
between the marshland and the water (represented in figure 2 as work gates 1, 2, 3, 
and 4). The conductivity of the interface is analogous to electrical conductivity. The 
conductivity of the interface is directly proportional to its length and to the thickness of 
the sheet of water flowing across the interface, as the conductivity of an electrical wire 
is proportional to its cross-sectional area. Thus, for a given hydraulic head, the rate 
of flow will be determined by the length of the land/water interface. In the case of 
work gate 6, the transport is impeded by a complex shoreline, as opposed to the first 
three gates. The more tortuous the path the flow must follow, the greater is the re- 
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sistance to flow. Thus, there is an inverse relationship between flow arid shoreline 
length or complexity for transport between the coastal waters and the bay-estuary sys- 
tem. 


Work gate 1 is bidirectional. Ebb tides remove material from the marsh end flood 
tides deposit material, with net transport being determined by the initial relative con- 
centrations of material. Work gate 3 is also bidirectional. A controlling factor in the 
transport of inorganic nutrients is the concentration of those nutrients in the bay water 
and in the interstitial waters of the marshland. Lower concentrations in the interstitial 
waters result in a fertilization of the marsh by the flood tide and little effect by the ebb, 
whereas a higher concentration of nutrients in the interstitial waters leads to removal 
of the nutrients by the tidal action. Work gates 2 and 4 are unidirectional, as rainwater 
falling on the marsh washes both detritus and inorganic nutrients across the shoreline 
into the bay. 

According to this reasoning, work gate 1, representing the interface conductivity, 
may be defined by the expression 


Wj = TSC i (1) 

where T is the mean tide range, S is the shoreline length, and Ci is a constant. 
The expression for work gate 2 is 


W 2 = RSC (2) 

where R is the amount of rainfall. The conductivity of the marsh-bay interface will 
therefore be defined by the composite function 

W mb = TSC i+RSC. = S (TC |+RC ) (3) 

The analogy of the ecosystem to an electronic circuit may be continued by com- 
paring the concentration or amount of detritus to voltage or electrical potential. The 
amount of detritus on the marshland is a function of the area of land interacting with 
the marine system and the primary production of the marshland. If there is more 
detritus in the bay than on the marshland subject to flooding, detritus will be left be- 
hind by the ebb tide. Conversely, higher levels of detritus on the marsh bottom will 
result in suspension of detritus by the flood tide and its removal to the bay subsystem 
by the ebb tide. Thus, the import of detritus, I (analogous to an electrical current), 
is related to the conductivity of the shoreline and detritus level, V, for the marsh and 
bay subsystems according to the simple equation 

I = ( Vir -V b )w mb (4) 

The factor V m is proportional to the amount of detrital material per unit area of marsh 
and the area subject to inundation or flushing by rainfall. It is normally significantly 
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greater than the "detritus potential" of the bay, V| j; although in some instances, 
may be greater. I positive indicates net flow from the marsh and I negative indicates 
net flow to the marsh. Assuming is negligible with respect to V m , equation (4) 
becomes 


1 = W mb V m (5) 

A similar relationship can be derived for export of detritus, E, from the bay to the 

coastal waters, where the detritus level, V c , is much lower than in the bays; the con- 

ductivity of the bay-coastal interface is defined as W^ c . 

E - w bc<W (6) 

The equilibrium detritus level, D, in the bay can then be written as 

D = f(I-E) (7) 

The function, f, includes consumption and sedimentation of detritus. 

The remaining factors of importance in figure 2 are the two switches controlled by 
salinity and temperature. Studies by the Louisiana Wildlife and Fisheries Commission 
(Barrett and Gillespie, 1973) indicate brown shrimp production is strongly affected by 
the number of hours water temperatures are below 20 # C after the first week of April, 
and that salinities over 10 ppt are required for a successful season, with 19 ppt close 
to the optimum salinity for the brown shrimp. To a first approximation, the switching 
functions might be represented as Gaussian curves 

K •• »2tt exp j^-b(p-p) “/£' J (8) 

where p represents the optimum salinity or temperature and $ defines the broadness of 
the curve, and hence the steepness of the switching function. An approximation of this 
type, although obviously crude, requires careful selection of the salinity and tempera- 
ture values to be used, as shrimp growth and production are related to these factors 
in a seasonal manner. 

If one makes the simplifying assumption that shrimp production is not affected sig- 
nificantly by predation and that there are no other factors important to the inshore pro- 
duction of shrimp, then it is possible to describe shrimp production, P, in Louisiana 
inshore waters mathematically as 

P - 8(D)K tGm pK sa j^ n £ t y (9) 

where g(D) is a presently undefined mathematical function. Because D is essentially 
determined by I (the flow of detritus from the marsh), shrimp production is closely 
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related to I, and hence determined by marsh productivity, shoreline complexity, and 
tidal and rain-induced flow. 

The commercial shrimp harvest is an indication, although probably not rerfect, of 
shrimp productivity. Because the productivity, P, is related to the import of detritus, 

I, a mathematical relationship should exist between the factors determining I and the 
commercial harvest. These relationships should be apparent as significant correlations 
and should make possible a predictive model. 


DATA 


The data analyzed fall into three categories: biological, physical, and geographi- 
cal . The biological data consist of the average inshore commercial shrimp catch for 
the years 1967 through 1972 reported by the National Marine Fisheries Service and 
tabulated as shrimp yield per acre in Barrett and Gillespie (1973). Pink, brown, and 
white penaeid shrimp contribute to these totals. Because of the intense fishing pressure 
and the economic factors involved, the yield data are very closely related to shrimp 
production and are used here as a measure of production. Temperature and salinity 
data reported oy Barrett (1971) were averaged for the period of April through August, 
1968. Mean tidal ranges listed for various points along the Louisiana coast in the 
National Ocean Survey Tide Tables were averaged for each geographic unit into which 
the coast was divided for this study, except for one area, where only a rough estimate 
of tidal range was available. The biological and physical data are presented in table I. 
The geographic data are derived from Landsat images of the Louisiana coastal region. 
The dates and scene identification codes of the Landsat data used in the study are listed 
in table II. 

The Louisiana coast has been divided into nine geographic units, shown in figure 3. 
The numbered units correspond to (1) Lakes Pontchartrain and Maurepas; (2) Lake 
Borgne and Chandeleur Sound; (3) Breton Sound; (4) the southern portion of the Missis- 
sippi River Delta; (5) Barataria Bay; (6) Terrebonne and Timbalier Bays; (7) the area 
extending from west of Terrebonne Bay to Atchafalaya Bay, including Caillou Bay; 

(8) Atchafalaya Bay through Vermilion Bay; and (9) from Vermilion Bay through Cal- 
casieu Lake. In general, the northern limit of the study area was taken to be the Intra- 
coastal Waterway. The unleveed marshes west of Lake Salvador and west and north of 
Lakes Pontchartrain and Maurepas were also included, whereas the leveed areas south 
of Lake Pontchartrain and along the Mississippi River were excluded. The nine areas 
include nearly all of the shrimp nursery grounds and inshore shrimping area of the state. 

The Landsat data available for use consist of photograph-like images and computer- 
compatible tape recordings of earth scenes. Each scene consists of registered images 
in four spectral bands: green, red, and two near infrared bands. Figure 4 is an example 
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TABLE I. - SHRIMP YIELD AND PHYSICAL FACTORS 


Geographic 

unit 

Shrimp 

yield a 

(kg/ha) 

Salinity* 3 

d>pt) 

c 

Temperature 
__ (°C) „ 

d 

Range 

(ft) 

1 

0.29 

4. 60 

28.00 

0.50 

2 

2.17 

13.14 

28. 01 

1.17 

3 

9.82 

15.14 

27.87 

1.33 

4 

7.88 

2. 12 

28.30 

1.25 

5 

35.57 

14.52 

27. 06 

1.03 

6 

45.16 

20.33 

27.38 

1.28 

7 

30.64 

15.84 

28.00 

1.65 

8 

1.10 

2.12 

28.30 

1.72 

9 

16.82 

9.20 

28. 65 

2.25 


a. From Barrett and Gillespie (1973) data on 1967-1972 commercial catch. 

b. From Barrett (1971) data for April through August, 1968, except value for unit 1 is 

from Stern and Atwell (1968) data for June and July, 1968, and value for unit 4 
is estimated to be the same as that for unit 8. 


c. From Barrett (1971) data for April through August, 1968. 


d. From National Ocean Survey Tide Tables. 



TABLE 11. - LANDSAT DATA USED IN ANALYSIS 


Date 

Scene identification 

Areas 

covered 

12 Oct. 1975 

2263 - 15491 

1, 2 

21 Oct. 1975 

5185 - 15323 

2, 3, 5 

21 Oct. 1975 

5185 - 15325 

3, 4, 5, 6 

25 Sept. 1975 

2246 - 15550 

o~i 

Ci 

OO 

25 Sept. 1975 

2246 - 15553 

6, 7, 8 

26 Sept. 1975 

2247 - 16005 

8, 9 

27 Sept. 1975 

2248 - 16063 

9 
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of the images obtainable from Landsat. Each image is composed of individual sample 
cells, referred to as picture elements. A picture element is approximately 57 meters 
wide (approximately the east-west direction) and 79 meters high (approximately the 
north-south direction). Computer analysis of the data permits conversion of the color 
data of the original scene into various thematic renditions. Using a standard computer- 
implemented image classification procedure, referred to as Water Search, land-water 
thematics were produced for the coastal region. A second computer program was used 
to geographically reference the data to the Universal Transverse Mercator System. 

This analysis permits the translation of points located on a map into the satellite co- 
ordinate system defining the thematic. The geographic coordinates defining the boun- 
daries of the nine geographic units were read from standard maps and translated into 
the satellite system to define the same boundaries in the thematics. A third computer 
program was then used to measure the shoreline length within the boundaries defining 
each geographic unit and to compute the shoreline density for each resolution element 
within each unit. A detailed description of the processing required to develop the shore- 
line length measurement is contained in Faller (1977). 

Shoreline density is defined as the length of shoreline per unit area. It is measured 
in the computer by scanning a window of predetermined size over the thematic (still in 
digital, computer-compatible form), accumulating the shoreline length within that win- 
dow, and then dividing the length by the area of the window. The shoreline density with- 
in the window is recorded for the reference picture element at the center of the window. 
The number of elements falling within each predefined range of densities is accumulated, 
and from this the total area described by each density range within a geographic unit is 
computed. The limits of the ranges of shoreline density used to analyze the Landsat 
data for this study are found in table III. The first range had less than 1. 0 km shore- 
line/km 2 , the second between 1. 0 and 2. 0 km, and so on. The window used was six 
picture elements high and eight wide, yielding ? nearly square window of about 465 
meters per side. 


TABLE III. - SHORELINE DENSITY RANGE 

Range 

1 2345678910 

Classification 

Minimum density — 12345679 11 

Maximum density 1 2345679 11 — 


Note: Density is in kilometers of shoreline per square 
kilometer of unit area. 
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Approximately 10 man-weeks of effort were expended in processing the Landsat 
data to generate the shoreline length and density products from the original data. 

The shoreline length and areas of land and water derived from the Landsat imagery 
are presented in table IV for each geographic unit. The normalized area of each unit 
falling within each shoreline density category, the normalized shoreline length, and 
the shoreline complexity factor are the data used in the study. The shoreline complex- 
ity factor is defined as the ratio of the actual shoreline length to that which would result 
from all the water being contained in a single circular lake, and is computed from the 
equation Q = hS/VTK, where S is the actual shoreline length and A is the area of water, 
both determined from the satellite data. It has a minimum value of 1.0, for the case 
of a perfectly circular lake, and increases as the number of lakes of decreasing size 
increases, or as the number of islands increases, or as the shoreline becomes con- 
voluted with small bays. Shoreline length is normalized by dividing the length measure- 


TABLE IV. - BASIC LANDSAT MEASUREMENTS 


Geographic 

Unit 

Land 

area 

(km 2 ) 

Water 

area 

(km 2 ) 

Shoreline 

length 

(km) 

1 

521 

494 

619 

2 

751 

4344 

2723 

3 

1080 

1512 

3419 

4 

225 

1044 

1318 

5 

2973 

2243 

8023 

6 

1268 

1419 

4005 

7 

1503 

1806 

4292 

8 

1125 

1538 

1642 

9 

2825 

2516 

4446 
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ment by the total area of the geographic unit. The shoreline density measurements 
were normalized by dividing the area classified into each density range by the active 
area of the geographic unit, defined as the total area of the unit from which is sub- 
tracted the area of land not falling within 230 meters of the shoreline. The area of 
land more than 230 meters from the shoreline was a byproduct of the shoreline density 
measureim nt, as the reference element at the center of the scanning window must be 
at least 230 meters from water for the window to be completely filled with land. These 
data are presented in table V. Also included in table V is the total area, having 
shoreline density greater than 5. 0 km shoreline/km 2 (range 6 and greater), normalized 
by tne active area. 


ANALYSIS 


The initial analytical effort was to examine correlations between the various param- 
eters described in previous sections and the shrimp productivity as indicated by the 
commercial shrimp yield. Linear correlation coefficients were computed according to 
the relation 

n 

I t (x r x) (y i -y'> 


where x and y are the mean values of the parameter being tested and the shrimp yield, 
and Sx and sy are the sample deviations for the parameter and the shrimp yield. The 
correlation coefficients are listed it; table VI, along with the significance level cf the 
correlations. Shrimp yield was plotted against some of the parameters and regression 
lines were computed from the data. Figure 5 is a graph of shrimp yield as a function 
of the normalized shoreline length, S; figure 6 shows yield as a function of shoreline 
complexity, Q; and figure 7 is a plot of normalized area having a shoreline density be- 
tween 5.0 and 6.0 km shoreline/km 2 . Statistical models 1, 2, and 3, found in table VII, 
are the least square error relationships between the shrimp yield and the respective 
parameters. Boot mean square (rms) deviations for these models are 8. 24, 8.62, and 
5. 73 kg/ha, respectively. The range of recorded shrimp yields is 0. 29 to 45. 16 kg/ha. 
Statistical model 4 relates the area falling into ranges 6 through 9 (>5. 0 km shoreline/ 
km 2 ) to shrimp yield. The rms deviation for this model was 6.49, not as good as that 
for model 3. 

The theoretical energy flow analysis discussed earlier suggests that shrimp pro- 
duction should be related to the product of the area producing detritus transported into 
the bay-estuary system, the tidal range, and the length of the shoreline (work gate 1, 
equation 1). It is assumed that the area producing detritus for export to the marine en- 
vironment is proportional to the area having a shoreline density greater than 5. 0 km 
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TABLE V. - DERIVED LANDS AT MEASUREMENTS 
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TABLE VI. - SHRIMP YIELD CORRELATION ANALYSIS 





Correlation 

Significance 

Factor 


Symbol 

coefficient 

level 

Normalized shoreline length 

S 

0.76 

0.98 

Shoreline complexity factor 

Q 

0.74 

0.98 

Land/water ratio 


L/W 

0.43 

<0.80 


r l 

a l 

0.26 

<0.80 


2 

a 2 

0.37 

<0.80 


3 

a 3 

0.23 

<0.80 

Normalized area for 

4 

a 4 

0.73 

0.98 

shoreline density i 

range ' 

5 

a 5 

0.78 

0.99 


6 

a 6 

0.83 

0. 995 


7 

a 7 

0.81 

0.99 


8 

a 8 

0.76 

0.98 


<9 

a 9 

0.28 

<0.80 

Tide 


T 

0.11 

<0.80 

Salinity 


0 

0.67 

0.95 

Temperature 


K 

-0.58 

0.90 

Normalized area for 
the sum of shoreline 
density ranges 6-9 


E 

0.81 

0.99 



TE 

0.77 

0.98 



SE 

0.78 

0.99 



STE 

0.84 

0.995 





1 



Figure 5. - Shrimp yield plotted against normalized 
shoreline length (S). 



Figure 6. - Shrimp yield plotted against shoreline 
complexity factor (Q = HS/v^ttA)* 
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NORMALIZED AREA 


Figure 7. - Shrimp yield plotted against normalized 
area with shoreline density greater than 5. 0 
and less than 6. 0 km shoreline/km 2 , ag. 
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TABLE VH. - SHRIMP YIELD MODELS 


RMS deviation 

Model Equation (kg/ha) 


1 

Y = 35.56S - 20.06 


8.24 

2 

Y = 1.08Q - 7.37 


8.62 

3 

Y = 818.8 a 6 - 12.49 


5.73 

4 

Y = 344 I - 7.13 


6.49 

5 

Y = 177.5 SET - 3.56 


4.94 

6 

Y = (186.7T - 11.35) SI - 

■ 3.64 

4.93 

7 

Y = (1. 849 x 10 4 T - 2. 026 x 10 3 ) SI - 208. 2 K a 

4.36 


where K 0 = 27>89 /5i exp 

|-ls[(a-39.57)/27.89] 2 } 



, 1 , N 

RMS deviation = %l — T 

y n 

! (Y i"^i )2 



Y = Predicted commercial shrimp harvest 
S = Normalized shoreline length 
Q = Shoreline complexity factor = ^SA/trA* 

A w = Area of water 

a 6 = Normalized area with shoreline density greater than 5. 0 and less than 
6. 0 km shoreline Am 2 

I = Normalized area with shoreline density greater than 5. 0 km shoreline/ 
km 2 

T = Mean tidal range 
o = Salinity 
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shoreline Am 2 . The product of these parameters was computed, and its correlation 
with shrimp yield was then determined to be 0. 84, thi highest of any of the parameters, 
and significant at a 99. 5-percent level. Model 5, the first statistical model based on 
ecological principles, was developed from this product and had an rms deviation from 
file actual shrimp yield of 4. 94 kgAa, a significant improvement over the first four 
models. 

The theoretical discussion also suggests that the product of shoreline length and the 
detritus-producing area should be related to shrimp production (rain-driven detritus 
export: work gate 2, equation 2). The correlation coefficient for this factor is 0. 78, 
but its inclusion improves the agreement between the prediction and the reported shrimp 
yield only slightly, resulting in an rms deviation of 4. 93 kgAa. 

As stated in the theorfc' 4 oal discussion, salinity and temperature are controlling 
factors in determining shrimp production. Models were generated which incorporated 
these two parameters in the form of Gaussian switches. These switches were repre- 
sented by factors 

K t = >2777 exp | ju-t)/f, t J 

K a = ex p j- 1 ; [('^)/«o]‘| 

where t, q, S t , and £ 0 are fitting parameters. The first two determine the point at 
which the switch is completely closed (i.e. , the mathematical value is maximum), 
whereas the latter two determine the steepness of the switching function. The tem- 
perature switch did not contribute to the model, and in fact worsened the agreement 
between predicted and measured shrimp yield values. The salinity switch did improve 
the agreement, with model 7 giving an rms deviation of 4. 36 kgAa. The model 7 pre- 
diction and reported shrimp yield are shown in figure 8. 


DISCUSSION 


The high correlations between the shoreline-related parameters (i.e., shoreline 
length and complexity factor and areas with high shoreline densities) and the shrimp 
yield are very convincing arguments in support of the hypotheses that the shoreline is a 
controlling factor in the flow of energy, stored in the form of detritus and its associ- 
ated microorganisms, from the marsh ecosystem into the bay-estuary ecosystem. The 
accuracy of model 5, based simply on the product of shoreline length, tidal range, and 
area with high shoreline density, supports the hypothesis that the detritus level in the 
bay-estuary subsystem is the main factor in determining the relative long-term shrimp 



SHRIMP YIELD KG/HA 


S.i 

K 



GEOGRAPHIC UNIT 


Figure 8. - Comparison of reported commercial shrimp 
landings with predicted yield from model 7 
based on shoreline length, area with shore- 
line density greater than 5.0 km shoreline/km 2 , 
tidal range, and salinity. 
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productivity along the Louisiana coast and that these levels arc controlled by the area 
producing detritus for export from the marsh and the interface between marsh and ma- 
rine subsystems. 

The small improvement of model 6 over model 5 (i. e. , the i nprovement resulting 
from the incorporation of the second work gate) indicates that the tide-independent flow 
of detritus from the marsh is not significantly different from the tide-dependent flow. 

The coefficients associated with the two types of flow (tide anti rain-driven) were very 
highly corr lated; in fact, the coefficient associated with the rain-driven flow is nega- 
tive, indicating that rain-driven flow of detritus is from the water to the land, an un- 
tenable conclusion. It is, of course, always dangerous to attempt to attach meaning to 
regression weights. This is particularly true when the two variables under considera- 
tion are highly correlated, as is the case with the two types of flow. The only accept- 
able conclusion is that, given the small tidal range typical of the study are-', he tide- 
dependent flow and the tide -independent flow are not statistically separable in terms of 
their effect on shrimp production. A single work gate would therefore suffice in place 
of work gates 1 and 2. 

The incorporation of the salinity switch resulted in some improvement in the pre- 
diction, although this improvement is small when one considers the importance of 
salinity in determining the success of one season as opposed to another. The effect of 
salinity is most apparent in geographic unit 6, the unit having the highest salinity and 
greatest shrimp yield. Model 6 predicted a yield of 39 kg/ha and model 7 predicted 
43 kg/ha, whereas the actual yield was 45 kg/ha. Models 6 and 7 differ only in the in- 
corporation Of the salinity switch. Selection of saPnity data for the analysis may in- 
fluence the significance implied by this analysis, as the year for which data were avail- 
able may not have been typical of the five years over which the yield data were averaged. 
If salinity data for all five years were available, the importance of salinity in deter- 
mining shrimp productivity might be more apparent in the model results. Another con- 
sideration is the salinity sampling locations. The points at which the measurements 
were made may not completely represent the nine geographic units, as some important 
portions of a given area may have a much different salinity from any of the points 
sampled in the unit. A : lathematicai form other than the Gaussian expression may also 
be more appropriate. 

The failure of the temperature switch in the model is probably due to the fact that 
the temperature data used represented the season average for a single year of the five 
for which the shrimp data were accumulated. The model would probably be improved 
by replacement of these data with the average number of hours water temperature was 
below 20 *C after the first week of April for each of the five seasons for each geographic 
unit. 


It is apparent that the predictions for geographic unit 3, Breton Sound, arc sig- 
nificantly higher than the actual yield. According to the model, tnis area should be very 
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productive, although in fact the shrimp yield is relatively low. A possible explanation 
of this involves migration patterns of the larval and postlarval shrimp offshore. The 
higher yield areas open directly on the Gulf of Mexico west of the Mississippi River, 
from the mouth of which flows tremendous volumes of fresh water, whereas Breton 
Sound is partially cut off from the open Gulf by the discharge of the river. It has been 
suggested by Barrett (personal communication) that the discharge of cold, fresh water 
by the Mississippi River may serve as a barrier interfering with the migration of post- 
larval shrimp found offshore, preventing them from entering the inshore waters of 
Breton Sound. 

Unfortunately, there are no adequate data to rigorously test the statistical signifi- 
cance of the models. Further research should be done to include in the analysis other 
coastal areas with large inshore and nearshore commercial shrimp harvests so that 
data points not used in the determination of the model coefficients can be used in testing 
the model. Despite the absence of independent test data, the accuracy of the models is 
such that the basic hypothesis (that the shoreline is a controlling factor in the production 
of shrimp through its regulation of the transport of detritus from the marshlands to the 
bays and estuaries, which constitute the shrimp nursery grounds and inshore harvest 
area) is strongly supported. 

A more rigorous analysis of the satellite data would probably result in slightly 
better agreement between yield predictions and actual yield, and would significantly 
improve the physical interpretation of the model. The current analysis assumed that 
the entire land area near the shoreline produced detritus uniformly. A more detailed 
analysis of the satellite data would differentiate vegetation species, using a currently 
available technique, and possibly estimate vegetation density (stems per square meter), 
using a teumique under development; the resulting thematic would be combined with a 
measure of distance from shore (a measurement technique that is now available) and 
provide a better estimate of detritus production. With the development of the mathe- 
matical relationship between detritus production and shrimp yield, the assignment of 
economic value to each unit of marsh land in terms of the shrimp industry would be 
possible, and the impact on shrimp production of a proposed modification of the marsh 
could be predicted. Trends of changes taking place in the marsh, whether natural or 
anthropogenic, could be analyzed in the light of this relationship to forecast possible 
changes in future shrimp production. Satellite data can be processed quickly at a 
reasonable cost to survey wide areas, even in remote coastal wetlands. The result is, 
quite possibly, a very powerful tool for resource management. 


CONCLUSIONS 


The flow of detritus from the marsh to the bays and estuaries of Louisiana appears 
to be a critical factor in determining inshore shrimp productivity. The commercial 


harvest of shrimp reported over a five-year period is highly correlated with shoreline 
length and complexify, and the area of land and water separated by a complex shore- 
line. Remote sensing techniques were used to develop a quantitative assessment of 
coastal shoreline features. Computer analysis of Landsat MSS data generated a map 
of the Louisiana coastal wetlands coordinated with shoreline length and density measure- 
ments. The techniques provided a current and accurate mapping of an area typified 
by constantly changing geography, at a very reasonable cost, demonstrating their po- 
tential for wide-area monitoring applications. 

The geographic data derived through these remote sensing techniques were used in 
correlation studies to examine the relationships between them and the commercial 
shrimp harvest. The geographic data were then used in several statistical models in 
conjunction with other physical data to predict the long-term average harvest. 

Landsat-based measurements of shoreline length and area of land and water having 
more than 5.0 km shoreline/km 2 were developed and used with published tidal ranges 
and salinities to predict the commercial shrimp yield for nine geographic units along 
the Louisiana coast; the rms deviation from the reported yield was 4.36 kg/ha over a 
range of 0. 29 to 45. 16 kg/ha. The mathematical model relating these parameters and 
the shrimp yield is consistent with an energy flow model describing the interaction of 
detritus -producing marshlands with shrimp nursery grounds and inshore shrimp fishing 
areas. The analysis of the geographic and physical parameters with the shrimp yield 
data thus supports the hypothesis that the shoreline is a controlling factor in the pro- 
duction of shrimp through its regulation of the transport of detritus. 

Day et al. (1973) observed that the most productive area of the estuary studied 
in the Barataria Bay region was along the marsh-water interface. They made the 
following observations: marsh grasses near the shore are often twice as high as those 
on the interior marshlands; the highest standing crops of marsh macrofauna and meio- 
fauna occur along the shore; standing crops of organic matter and meiobenthos in the 
submerged sediments are higher near shore; and benthic populations are densest near 
the shore. Day et al. state: "These factors suggest that overall marsh production 
will increase as the amount of marsh edge habitat is increased. The familiar picture 
of salt marshes with many twisting and dendritic channels probably reflects a tendency 
of the estuary system to develop maximum production." What Day et al. observed in 
a broad range of species of flora and fauna over a very restricted area, this study has 
demonstrated quantitatively for a single organism (the penaeid shrimp) of bays and 
estuaries along the entire Louisiana coast. 


National Space Technology Laboratories 

National Aeronautics and Space Administration 

N3TL Station, Mississippi 39529 January 17, 1979 


25 



REFERENCES 


Barrett, Barney B. : Cooperative Gulf of Mexico Estuarine Inventory and Study, Louisi- 
ana. Phase II, Hydrology, and Phase III, Sedimentology. La. Wildl. and Fish. 
Comm., 1971. 

Barrett, Barney B. , and Marilyn C. Gillespie: Primary Factors which Influence Com- 
mercial Shrimp Production in Coastal Louisiana. Technical Bulletin No. 9. La. 
Wildl. and Fish. Comm. , 1973. 

Carter, M. R. , L. A. Burns, T. R. Cavinder, K. R. Dugger, P. L. Fore, D. B. 
Hicks, H. L. Revells, and T. W. Schmidt: Ecosystems Analysis of the Big Cypress 
Swamp and Estuaries. U.S. Environ. Prot. Agency Region IV, 1973. 

Day, J. W. , Jr., W. G. Smith, and C. S. Hopkinson, Jr.: Some Trophic Relation- 
ships of Marsh and Estuarine Areas. Proc. Coastal Marsh and Estuary Mgt. 

Symp. La. State Univ. , pp. 115-135, 1972. 

Day, J. W. , W. G. Smith, P. R. Wagner, and W. C. Stowe: Community Structure 
and Carbon Budget of a Salt Marsh and Shallow Bay Estuarine System in Louisiana. 
Publication No. LSU-SG-72-04. La. State Univ. , 1973. 

Faller, K. H. : A Procedure for Detection and Measurement of Interfaces in Remotely 
Acquired Data Using a Digital Computer. NASA TR R-472, 1977. 

National Ocean Survey: Tide Tables. High and Low Water Predictions, East Coast of 
North and South America Including Greenland. U.S. Dept. Comm. , Nat. Oceanic 
and Atmos. Admin., 1978. 

Odum, W. E., J. C. Zieman, and E. J. Heald; The Importance of Vascular Plant 
Detritus to Estuaries, pp. 91-114. Proc. Coastal Marsh and Estuary Mgt. Symp., 

R. H. Chabreck, ed. La. State Univ. , 1972. 

Palmisano, A. W. , Jr.: Plant Community-Soil Relationships in Louisiana Coastal 
Marshes. PhD Dissertation, La. State Univ. , 1970. 

Schelske, C. L. and E. P. Odum: Mechanisms Maintaining High Productivity in 
Georgia Estuaries. Proc. Gulf Carib. Fish. Inst. , 1975-80. 1961. 



Stem, D. H. , and B. H. Atwell: A Summer Limnology Study of Lake Pontchar train, 
Louisiana. La. Water Resources Research Inst. , Tech. Rpt. No. 3, La. State 
Univ. , 1968. 

Teal, J. M.: Energy Flow in the Salt Marsh Ecosystem of Georgia. Ecology, 43(4): 
614-624, 1962. 


27 


